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We construct inhomogeneous asymptotically AdS black hole solutions corresponding to the spon- 
taneous breaking of translational invariance and the formation of striped order in the boundary field 
theory. We find that the system undergoes a weak first order phase transition in both the fixed 
density and fixed chemical potential ensembles. The bulk solutions have striking geometrical fea- 
tures related to a magnetoelectric effect associated with the existence of a near horizon topological 
insulator. At low temperatures the horizon becomes highly inhomogeneous and tends to pinch off. 
Entropically and energetically, these solutions cannot serve as the end-state of the dynamical insta- 
bility, which strongly suggests that a transition to a solution localized in the transverse direction is 
taking place, in a process that potentially violates the cosmic censorship conjecture. 



I. INTRODUCTION AND SUMMARY 

The gauge-gravity duality describes phenomena in 
strongly coupled field theories via their relation to clas- 
sical or semi-classical gravitational systems. From the 
perspective of the boundary quantum field theory, this 
relation can be used to construct and study models for 
ill-understood phenomena which arise in such strongly 
coupled systems. On the other hand, the relation to local 
quantum field theory helps motivate and interpret new 
results in classical and quantum gravity. In this Letter we 
apply holography to study the spontaneous breaking of 
translation invariance and the formation of striped order. 

Stripes are known to form in a variety of strongly cou- 
pled systems, from large N QCD [1, 2] to systems of 
strongly correlated electrons (for a review see [3]). The 
formation of stripes and the associated reduced dimen- 
sionality are speculated to be related to the mechanism 
of superconductivity in the Cuprates [4]. It is therefore 
useful to study striped phases in the holographic context. 

Besides its interest in the boundary theory, this study 
has an intrinsic interest in the bulk gravitational con- 
text. Our model describes a black hole whose instability 
to the formation of inhomogeneous structures resembles 
the black string instability [5] . That instability is known 
to lead to horizon pinch off and the violation of cosmic 
censorship in five dimensions [6] and above [7]. Similar 
to the black string case, the inhomogeneous solutions we 
construct are found to have smaller entropy than the ho- 
mogeneous solutions and a larger mass than the critical 
black hole, and thus cannot serve as the end-state of the 
instability in four dimensions. This suggests that the dy- 
namical formation of a naked singularity is possible in 
the decay process, and thus may provide an example of 
cosmic censorship violation in four dimensions. 

Our study is facilitated by a numerical solution of 
the set of coupled non-linear Einstein and matter equa- 
tions in the bulk, which exhibit a normalizable inhomo- 
geneous mode. Previous studies of inhomogeneous so- 
lutions in asymptotically AdS spacetimes concentrated 
on non-normalizable modes (i.e. explicit rather than 



spontaneous breaking of translation invariance) or the 
study of co-homogeneity one solutions [8-11], where one 
of the translational Killing vectors is replaced by a helical 
Killing vector. More recently, a probe model was shown 
to have a magnetic field induced lattice ground state [12]. 
In contrast to the above, our solutions backreact on the 
geometry, are co-homogeneity two, and exhibit sponta- 
neous breaking of translation invariance below a critical 
temperature. 

We describe striking bulk and boundary properties of 
our bulk solutions, including frame dragging effects, the 
magnetic field, the curvature and the geometry. Some 
of the features can be understood as the emergence of 
a near horizon region which acts as a bulk topological 
insulator. The magnetoelectric effect is then responsible 
for the patterns we observe for the bulk magnetic field 
and vorticity. 

These features are analyzed as a function of tempera- 
ture. In particular, we find that the horizon of the black 
hole develops a "neck" and a "bulge" in the transverse 
direction which shrink with temperature, such that the 
ratio of their sizes contracts as fast as <~ T 1 / 2 . Simultane- 
ously, the proper length of the horizon in the transverse 
direction grows at a rate ~ 1/T 01 . However, the curva- 
ture remains finite, and its maximal value, occurring at 
the bulge, tends to a constant in the limit T — > 0. 

The bulk black hole solutions give rise to the holo- 
graphic stripes on the boundary, characterized by non 
vanishing momentum and electric current and modula- 
tions in charge and mass density. Starting small near T c , 
the amplitudes of the modulations grow steadily at lower 
temperatures, approaching finite values at T — > 0. 

We study the thermodynamics of the system by con- 
structing phase diagrams in various ensembles. We dis- 
cover a weak first order phase transition both in the 
canonical (fixed charge) and in the grand canonical (fixed 
chemical potential) ensembles, to the accuracy of our 
numerics. On the other hand, analysis in the micro- 
canonical ensemble, most relevant to dynamical evolu- 
tions of isolated gravitational systems, shows that the 
smooth inhomogeneous solutions we find cannot be the 
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end-state of the instability. This signals that the dynam- 
ical system proceeds violently to a different state, likely 
with the emergence of naked singularity in the process. 



II. THE HOLOGRAPHIC SETUP 

The Lagrangian describing our coupled system is [13] 

c = \r - \d^d^ - \f^f^ - - d nt , 
V(tp) = -6 + ^m 2 ip 2 , 

where R is the Ricci scalar, F^ v is the Faraday tensor, 
Hint describes the axion coupling and g is the determi- 
nant of the metric. We use units in which the AdS radius 
I 2 = 1/2, Newton's constant 8ttGn = 1, and c = % = 1, 
and choose m 2 = —2 and c\ = 4.5. Pcrturbative in- 
stabilities towards the formation of charge and current 
density waves were identified in [13] for a range of wave 
numbers and temperatures 1 . We note the appearance of 
axion electrodynamics in the bulk theory. It is curious 
that inhomogeneous instabilities (see also [14]) seem to 
involve the topology of the bulk fields in an essential way, 
though the analysis performed to discover the instability 
is local in nature. We leave this mystery for future work. 

In this Letter we investigate the end-point of the in- 
stability driven by the wave number k c ~ 0.75 with the 
largest critical temperature T c ~ 0.012 [13]. This state 
is a co-homogeneity two solution, thus we construct the 
family of stationary solutions that emerge from the crit- 
ical point, assuming all the fields to be functions of the 
radial coordinate r and one spatial coordinate x. 

Our ansatz includes the scalar field ip{r,x), the gauge 
field components A t (r, x) and A y (r, x) and the metric 

ds 2 = -2r 2 f(r)e 2A ^dt 2 + 2r 2 e 2C ^ (dy - W(r, x)dtf 

+ e2Bir - X) (2^ 2 W) +2r2dx2 )> (2) 

where for the sake of convenience we included in the def- 
inition of the metric functions the factor /(r) character- 
izing the metric of the AdS Reissner-Nordstrom (RN for 
short) solution, with horizon at r = r : 

'<'>-'-(' + + 

The inhomogeneous solutions reduce to the RN solution 
above the critical temperature. 



1 An interesting application of the instability in this model has 
appeared very recently [15]. 



The conformal in r, x plane ansatz (2) is convenient 
in constructing co-homogeneity two solutions. With this 
ansatz, the Einstein and matter equations reduce to seven 
coupled elliptic equations and two constraint equations. 
Moreover, the constraint system can be solved elegantly 
using its similarity to a Cauchy-Riemann problem [16]. 

The boundary conditions we impose correspond to 
regularity conditions at the horizon and asymptotically 
AdS conditions at the conformal boundary. With these 
boundary conditions, the set of solutions we find depend 
on three parameters: the temperature T, the chemical 
potential p and the periodicity in the x direction L. Us- 
ing the conformal symmetry inherent in asymptotically 
AdS spaces, the moduli space of solution depends only on 
the two dimensionless combinations of these parameters. 
Furthermore, here we discuss only the dominant critical 
mode that becomes unstable at the largest temperature 
T c , so that L = 2n/k c . 

On the spatial boundaries it is useful to impose "stag- 
gered" periodicity conditions. Using two reflection sym- 
metries which are preserved by the form of the unstable 
perturbation, one can reduce the numerical domain to a 
quarter period and impose 

d x <P(x = 0) - 0, V(z = \ ) = o, 

h(x = 0) = 0, d x h(x = -) = 0, 

d x g(x = 0) = 0, d x g(x=j)=0, 

where h represents the fields A y and IT, and g refers 
collectively to A, B, C and A t . 

The elliptic equations derived from (1) arc discretized 
using finite difference methods and are solved numeri- 
cally by a straightforward relaxation with the specified 
boundary conditions. In this method the equations are 
iterated starting with an initial guess for all fields, un- 
til successive changes in the functions drop below the 
desired tolerance. We verify that the remaining two con- 
straints are satisfied by those solutions. Full details of 
our procedure are given in the upcoming [17]. 

III. THE SOLUTIONS 

A convenient way to parametrize our inhomogeneous 
solutions is by the dimensionless temperature 9 = T/T c , 
relative to the critical temperature T c . Our method al- 
lows us to find solutions in the range 0.003 < 9 < 0.9. 

Bulk Geometry. For subcritical temperatures, as 
we descend into the inhomogeneous regime, the metric 
and the matter fields start developing increasing vari- 
ation in x. Fig. 1 displays the metric functions for 
9 ~ 0.11, over a full period in the x direction. The matter 
fields have qualitatively similar behaviour. The variation 
of all fiel ds is ma ximal near the horizon of the black hole 
at p = \Jr 2 — r 2 = 0, and it gradually decreases toward 
the conformal boundary, p — > oo. 
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FIG. 1: Metric functions for 9 ~ 0.11. Note that the metric 

functions A, B and C have half the period of W. The 
variation is maximal near the horizon, located at p = 0, and 
it decays as the conformal boundary is approached, when 
p — > oo. The matter fields (not shown) behave in a 
qualitatively similar manner. 



Many of the special features of the solutions we find 
are related to the presence of axion electrodynamics, the 
effective description of the electromagnetic response of a 
topological insulator, in the gravity action. In the bro- 
ken phase we have an axion gradient in the near horizon 
geometry, which therefore realizes a topological insulator 
interface 2 . The presence and the pattern of a near hori- 
zon magnetic field, summarized in the field A y , can be 
related to the magnetoelectric effect in such interfaces. 

In curved space the magnetic field is accompanied by 
vorticity, which is manifested by the function W. This 
causes frame dragging effects in the y direction. Test 
particles will be pushed along y with speeds W(r,x), in 
particular the direction of the flow reverses every half 
the period along x. The drag vanishes at the horizon 
and at the location of the nodes of W where x = Ln/2, 
for integer n (see Fig. 1). In general, the dragging effect 
remains bounded, the vector d t is everywhere timelike, 
and no ergoregion forms. 

The Ricci scalar of the RN solution is Rrn = —24, 
constant in r and independent of the parameters of the 
black hole. This is no longer true for the inhomoge- 
neous phases, where the Ricci scalar becomes position 
dependent. The right panel of Fig. 2 illustrates the spa- 
tial variation of the Ricci scalar, relative to the Rrn for 
9 ~ 0.003. The maximal curvature is always along the 
horizon at x — nL/2 for integer n. It grows when the 



2 It would be interesting to discuss localized matter excitations on 
the interface, especially fermions, along the lines of [18]. 
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FIG. 2: Left panel: The variation along x of the size of the 
horizon in the y direction includes alternating "necks" and 
"bulges". Right panel: Ricci scalar relative to that of RN 
black hole, R/Rrm — 1 for 8 ~ 0.003 over half the period. 
The scalar curvature is maximal along the horizon at the 
bulge x = nL/2 for integer n. 



temperature decreases and approaches the finite value of 
R ~ —94 in the small temperature limit. 

The left panel in Fig. 2 shows the variation of trans- 
verse extent of the horizon in the y direction, r y (x) = 
r exp[C(r , x)], along x for 9 ~ 0.003. Typically 
there is a "neck" occurring at x = nL/2 and a "bulge" 
at x = (n + 1) L/2, for integer n. Note that Ricci scalar 
curvature is maximal at the bulge and not at the neck as 
would happen, for instance, in the spherically symmetric 
black string case. The size of both the neck and the bulge 
monotonically decrease with temperature, however, the 
neck is shrinking faster. We find that the ratio scales as 
r ncck^ r buigc ^ qi/2 near ^he i ower en( \ f the range of 6's 

that we investigated. 

Another aspect of the geometry is the proper size 
of the stripe in the x direction at fixed r, l x {r) = 
Jq exp[£?(r, x)] dx. The proper length tends to the co- 
ordinate length as l/r 3 asymptotically as r — > oo, but it 
exceeds that as the horizon is approached. Namely, the 
inhomogeneous phase "pushes space" around it along x, 
resembling the "Archimedes effect". The proper length 
of the horizon is maximal and it grows as the tempera- 
ture decreases. We find that at small 9 the proper length 
of the horizon diverges as ~ 6*~ 01 . 

Boundary Observables. Near the conformal 
boundary the fields decay to their AdS values, and the 
subleading terms in their variation are used to define 
the asymptotic charge densities of our solutions. The 
subleading fall-offs of the metric functions in our ansatz 
determine the boundary stress-energy tensor, whereas 
the fall-offs of the gauge field determine the charge 
and current densities of the boundary theory. Finally, 
the subleading term of the scalar field near infinity 
determines the scalar condensate. 
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FIG. 3: Difference in the thermodynamic potentials between 

the inhomogeneous phase and the RN solution, plotted 
against the temperature. In both ensembles the differences 
are small, signalling a weak first order transition where the 
inhomogeneous solution dominates at small temperatures. 



For our inhomogeneous solutions we find that all 
charge and current densities are spatially modulated, ex- 
cept for (T xx ), which is constant, consistent with the con- 
servation of boundary energy-momentum. We define the 
total charges of a single stripe by integrating the charge 
densities over the full period L. These integrated quan- 
tities are charge densities per unit length in the transla- 
tionally invariant direction y. 



IV. THERMODYNAMICS 

We demonstrated that below the critical temperature 
T c there exists a new branch of solutions which are spa- 
tially inhomogeneous. The question of which solution 
dominates the thermodynamics depends on the ensem- 
ble used. We start our discussion by fixing the temper- 
ature and the total charge. This describes the physical 
situation in which the system is immersed in a heat bath 
consisting of uncharged particles. 

In the upper panel of Fig. 3 we plot the difference of 
the normalized total free energy, F = M — TS, between 
the two classical solutions as function of the tempera- 
ture T. In our ensemble the total charge N is fixed, and 
we use the scaling symmetry of the boundary theory to 
set N — 1, or in other words measure all quantities in 
terms of TV. As a result the free energy is a function of 
one parameter, the temperature T. The figure displays 
a first order phase transition, where the inhomogeneous 
solution dominates the thermodynamics below a critical 
temperature T c , below the temperature at which inho- 
mogeneities first develop. 

We note that the free energy above the critical tem- 
perature T c is almost degenerate with that of the homo- 
geneous solution, which results in uncertainty for deter- 
mining the value of the critical temperature. While our 
numerics is accurate to the level of under a percent, for 
the current parameter-setting these small inaccuracies re- 



FIG. 4: The entropy difference between the new 
inhomogeneous solutions and the RN background. We see 

that the difference is always negative. The mass of the 
solutions is always higher than that of the critical RN black 

hole. 



suit in large variations in the value of T c . However, we 
can state with high confidence that the phase transition 
exists, i.e that T c is finite and positive. 

If we fix the chemical potential instead of the charge, 
we discuss a situation where the system is immersed 
in a plasma made of charged particles. To study the 
thermodynamics we use the grand canonical free energy 
fi = M — TS — fiN, displayed in the lower panel of Fig. 
3. In this ensemble it is convenient to measure all quan- 
tities in units of the fixed chemical potential /i. Then, 
again, the free energy is a function of only the tempera- 
ture T. In the fixed chemical potential ensemble we find 
a similar first order transition, where the inhomogeneous 
charge distribution starts dominating the thermodynam- 
ics below the temperature where the inhomogeneous in- 
stability develops. 

The physical situation relevant to the study of the real 
time dynamics of the instability corresponds to fixing the 
mass and the charge. This is the microcanonical ensem- 
ble, describing an isolated system in which all conserved 
quantities are fixed. In this ensemble it is convenient to 
measure all quantities in terms of the (fixed) charge, and 
the remaining control parameter is then the mass M. We 
find that the mass of the inhomogeneous solution always 
exceeds that of the critical solution at T c . Furthermore, 
in Fig. 4 we display the entropy difference between the 
two classical solutions. The entropy of our inhomoge- 
neous solutions is always smaller than that of the RN 
black hole of the same mass. Thus a dynamical decay 
to our inhomogeneous solution is forbidden on both en- 
ergetic and entropic grounds. We believe that this may 
be a first indication that this model exhibits a violation 
of cosmic censorship in four dimensions. 
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